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9^. ABSTRACT 
<r ■ 

' We have built an infrared beam combiner for the GI2T/REGAIN interferometer of the Observatoire de la Cote d'Azur. 
CO ' The beam combiner allows us to record spectrally dispersed Michelson interference fringes in the near-infrared J-, H- 
' ' or K-bands. The beam combiner has the advantage that Michelson interferograms can simultaneously be recorded 
I ' in about 128 different spectral channels. The tilt of the spectrally dispersed fringes is a measure of the instantaneous 
J> optical path difference. We present the optical design of the beam combiner and GI2T/REGAIN observations of 
the Mira star R Gas with this beam combiner in the spectral range of 2.00 /im - 2.18 /im (observations on 22 and 
OO ' 25 August 1999; variability phase 0.08; V-magnitude approximately 6; seven baselines between 12 m and 24 m; 

reference stars Vega and (3 Peg) . The spectrograph of the beam combiner consists of an anamorphotic cylindrical 
lens system, an image plane slit, and a grism. The detector is a 256x256 pixel Rockwell PIGNIG array camera. The 
' shortest possible exposure time is 10 ms. A system of digital signal processors calculates the ensemble average power 
spectrum of the spectrally dispersed Michelson interferograms, the instantaneous optical path difference error, and 
' several other useful parameters in real time. From the observed R Gas visibilities at baselines 12.0 m, 13.8 m and 
13.9 m, a 2.1 fim uniform-disk diameter of 25.3 mas ± 3.3 mas was derived. The unusually high visibility values at 
' baselines > 16 m show that the stellar surface of R Gas is more complex than previously assumed. The visibility 
^ values at baselines > 16 m can be explained by high-contrast surface structure on the stellar surface of R Gas or 
other types of unexpected center-to-limb variations. The R Gas observations were compared with theoretical Mira 
star models^^^. We obtained the following results for R Gas at variablity phase 0.08: 
Angular Rosseland radius i?": 12.1 mas ± 1.7 mas 
•'H . Linear Rosseland radius R: 276 R0 ±66 Rq (using the HIPPARGOS parallax of 9.36mas±1.10mas) 
rN ' Effective temperature: 2685 K ± 238 K (derived from the angular Rosseland radius and the JHKLM-photometry) . 
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1. INTRODUCTION 

The GI2T/REGAIN interferometer^'^ is an optical interferometer of two 1.5 m telescopes. For this interferometer we 
have built an infrared beam combiner and used it to record spectrally dispersed infrared Michelson interferograms of 
the Mira star R Gas. In this paper we present the design of our new infrared beam combiner and R Gas observations 
obtained in August 1999. This beam combiner allows us to simultaneously record Michelson fringes in about 128 
different spectral channels from 2.00 to 2.18 /im. The tilt of the spectrally dispersed fringes is a measure of the 
instantaneous optical path difference. The R Gas observations reported in this paper are the first NIR interferometry 
observations with large 1.5 m telescopes. The aim of our Mira star project is to resolve the stellar disk of Mira stars, 
to reveal photospheric asymmetries and surface structures, and to study the wavelength and phase dependence of the 
diameter. Previous speckle or long-baseline interferometry observations of Mira stars were, for example, reported in 
Refs. 5-16. Theoretical studies^'^'^^^^^ show that interferometric diameter measurements can considerably improve 
our understanding of cool stellar atmospheres. 



2. INSTRUMENT 



Figure 1 shows the optical layout of our new NIR beam combiner built at the MPI for Radioastronomy. It consists of 
an anamorphotic lens system, an image plane slit, a grism, a shutter, a PICNIC detector, a data storage computer, 
and an array of digital signal processors for data processing. The DSPs can calculate the power spectrum and the 
instantaneous optical path difference in real-time. The typical frame rate is 3 frames per second (256x256 pixels per 
frame) . 
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Figure 1. Optical layout of the new NIR beam combiner (1 beam combiner lens with 750 mm focal length; 2, 
4 and 6 cylindrical lenses; 3 image plane slit and slit viewer camera; 5 grism; 7 PICNIC detector; anamorphotic 
compression factor 9). 



3. R CAS OBSERVATIONS AND DATA PROCESSING 

The first infrared observations with our new beam combiner mounted to the GI2T/REGAIN interferometer were 
carried out on August 22 and 25, 1999. The observational parameters of the R Cas and reference star observations 
are as follows: 



Projected R Cas baselines: 
Date of R Cas observations: 

Reference stars: 

Number of frames per baseline: 

Frame rate: 

Filter: 

Recorded wavelength band: 
Number of spectroscopic channels: 
Anamorphotic projection factor: 
Exposure time per frame: 
K-band seeing: 



12.01 m, 13.80 m, 13.88 m, 18.14 m, 19.79 m, 21.99 m, 23.79 m 
22 Aug. 1999 (baselines 12.01 m, 13.80 m), 

25 Aug. 1999 (baselines 13.88 m, 18.14 m, 19.79 m, 21.99 m, 23.79 m) 

Vega, (3 Peg 
500 

3 frames /s 

center wavelength 2.1 /im, FWHM bandwidth 0.3 /xm 

2.00 ^m - 2.18 ^m 

128 

9 

100 ms 
2 arcsec 



Figure || shows two of several thousand recorded interferograms of R Cas. The interferograms of R Cas and the 
reference stars were flatfielded and sky-subtracted, and the visibilities of each spectral channel and baseline were 
derived by calculating the ensemble average power spectra. The method described by Mourard et al.^° was applied 
to derive visibility ratios of the object and reference star observations and to obtain calibrated R Cas visibilities. 
Figure || shows the derived R Cas visibilities. The solid line is the visibility function of a uniform-disk of a star with 
25.3 mas diameter. The dashed curve is the uniform-disk visibility fit curve which is required for our observations 
with large 1.5 m telescopes. This dashed visibility fit function has no zeros since an average of many baselines is 
measured rather than only one single baseline if the telescope pupils are large. If the projected baseline from mirror 
center to mirror center is for example 10 m, all baselines between 8.5 m and 11.5 m are simultaneously measured. The 
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Figure 2. Two of several thousand recorded interferograms of the Mira star R Cas (spectral range 2.00 to 2.18 
/Ltm; baseline 13.80 m). In the left interferogram the fringes are nearly parallel to the dispersion direction since the 
optical path difference error is almost zero. In the right interferogram the fringes are tilted due to a small optical 
path difference error. 
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Figure 3. GI2T/REGAIN visibility measurements of R Cas with our new NIR beam combiner (baselines 12.01 m, 
13.80 m, 13.88 m, 18.14 m, 19.79 m, 21.99 m, 23.79 m). 



R Cas visibilities in Fig. |^ show that the center-to-limb variation of R Cas is not a uniform-disk function. The large 
visibility values at spatial frequencies larger than ~ 40 cycles per arcsec can be explained by high-contrast surface 
structure on the stellar surface of R Cas or by other types of unexpected center-to-limb variations. A uniform- 
disk fit diameter of 25.3 mas ± 3.3 mas is obtained if only the first three low-frequency visibility points (projected 
baselines 12.01 m, 13.80 m, and 13.88 m) are used for the uniform-disk fit. The error bars of our first IR observations 



are relatively large because the weather situation was not favorable and the control system (Mourard et al., this 
proceedings^^) of the GI2T/REGAIN (Mourard et al., this proceedings^^) was not fully operational. Under normal 
observing conditions the error bars will be much smaller. 

4. COMPARISON OF THE OBSERVATIONS WITH MIRA STAR MODELS 

4.1. The stellar filter radius and Rosseland radius of R Cas 

In this section we derive various types of angular diameters from the measured visibilities of R Cas by applying the 
different theoretical center-to-limb intensity variations (hereafter CLV) of different Mira star models (Bessel, Scholz, 
and Wood 1996 = BSW96\ Hofmann, Scholz and Wood 1998 = HSWgS^). From these angular diameters (Fig. |) 
and the bolometric flux, we derive effective temperatures (Fig. ||). All Mira star models used in this paper are from 
BSW96 (D and E series) and from HSW98 (P, M and O series). They were developed as possible representations of 
the prototype Mira variable o Ceti, and hence have periods P very close to the 332 day period of this star. They differ 
in pulsation mode, assumed mass M and assumed luminosity L. The BSW96 models differ from the more advanced 
HSW98 models with respect to the pulsation modelling technique. The five models represent stars pulsating in 
the fundamental mode (/; D, P and M models) or in the first-overtone mode (o; E and O models). Table H lists 
the properties of these Mira model series. Table || provides the link between the 22 abscissa values (model-phase 
combinations m) in Figs. § - ^, and the models. R-p denotes the Rosseland radius (below) of the non-pulsating parent 
star of the Mira variable (BSW96, HSW98). Table ^ additionally lists the variability phase, the relative Rosseland 
and stellar K-band radius, and the effective temperature. We compare predictions of these models for different phases 
and cycles with our observations. 

Monochromatic radius Rx and Rosseland radius R. We use the conventional stellar radius definition where the 
monochromatic radius Rx of a star at wavelength A is given by the distance from the star's center at which the optical 
depth equals unity {tx = 1). In analogy, the photospheric stellar radius R (Rosseland radius) is given by the distance 
from the star's center at which the Rosseland optical depth equals unity (tross !)• This radius has the advantage 
of agreeing well (see Table 6 in HSW98 for deviations sometimes occurring in very cool stars) with measurable 
near-infrared continuum radii and with the standard boundary radius of pulsation models with Toff oc {L/R'^y^^. 

Stellar filter radius R[. For the K-band filter used for the observations, we have calculated the theoretical CLVs 
of the above mentioned five Mira star models at different phases and cycles. The stellar radius for filter transmission 
f^ is the intensity- and filter- weighted radius Rf — J Rx I\ i\d\/ J Ix f\ rfA, which we call stellar filter radius R{ after 
Scholz & Takeda's definition^'^. In this equation Rx denotes the above monochromatic ta = 1 radius, Ix the central 
intensity spectrum and f^ the transmission of the filter. 

4.2. Angular R Cas radii 

The observed angular stellar K-band radii R^ of R Cas corresponding to the various model-phase combinations m, 
were derived by least-squares fits between the first three visibility values (baselines 12.01 m, 13.80 m and 13.88 m) 
shown in Fig. y and the visibilities of the different theoretical K-band CLVs of different model-phase combinations 
m. Additionally, the observed angular Rosseland radii R^ of R Cas corresponding to the model-phase combination 
m were derived from the observed angular stellar K-band radii R^ and the theoretical radius ratios Rm/ RK,m from 
Table || (Table || provides theoretical R and Rk values for each model-phase combination m) . The derived angular 
stellar filter radii i?K,in and angular Rosseland radii i?m are shown in Fig. |[ 



Table 1. Properties of Mira model series"'^'^ (see text) 
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R Gas: angular stellar K-band filter radii 
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R Gas: angular Rosseland radii 
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Figure 4. Observed angular stellar K-band radii (top) and Rosseland radii Rf^ (bottom) of R Gas derived 

from the measured R Gas visibilities by fitting the theoretical CLVs of the 22 model-phase combinations m. Table ^ 
shows the link between the abscissa values and the models and the variability phases. 



Table 2. Link between the 22 model-phase combinations m in Figs. ^ - |6|, and the models. The variability phase 
^vis, the Rosseland radius R and the K-band radius Rk in units of the non-pulsating parent star radius Rp (see 
BSW96), and the effective temperature Tcb{R) associated to the Rosseland radius are additionally given. 
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4.3. Linear R Cas radii 

We have derived linear stellar radii of R Cas (in units of solar radii) from the measured angular stellar K-band radii 
and Rosseland radii R^-^ (Fig. |) by using the R Cas HIPPARCOS parallax of 9.36 mas ±1.10 mas^^. Fig. | 
shows the obtained R Cas radii for all model-phase combinations m. The theoretical Rosseland radii of the M and 
of the P model series at nearly all available near-maximum phases (1.0, 3.0, 4.0) are close (within the error bars) to 
the measured Rosseland radii. The theoretical Rosseland radii of the first-overtone models E and O are clearly too 
large compared with measured Rosseland radii. The theoretical radii of the D model series are slightly too small 
compared with the measured Rosseland radii. The same conclusions are also valid for the linear stellar filter radii 
Rk (Fig. |). 

If we calculate average linear radii by averaging the radii obtained with all available near-maximum M models 
(i.e., m=17, 19) and/or near-maximum P models (i.e., m = 9, 11, 13, 15) we obtain: 

Average theoretical M model Rosseland radius: 315 R© 

Average measured M model Rosseland radius: 287R0 ±51 R© (obtained with m= 17 and 19) 

Average theoretical P model Rosseland radius: 263 R© 

Average measured P model Rosseland radius: 264Ro ±82 R© (obtained with m= 9, 11, 13 and 15) 



Average theoretical model Rosseland radius: 289 R© 

Average measured model Rosseland radius: 276 R© ±66 R© (m=:9, 11, 13, 15, 17, 19; M and P models) 



R Gas: linear stellar K-band filter radii 
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R Gas: linear Rosseland radii 
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Figure 5. Comparison of observed R Cas radii (crosses; see Sect. 4.3) and theoretical model radii (filled squares; 
D, E, P, M and O model for different phases; see Table |) plotted versus all 22 model-phase combinations m (top: 
stellar K-band radii i?K m! bottom: Rosseland radii i?m)- 



R Gas: effective temperatures 
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Figure 6. Comparison of model effective temperatures Toff^m (plotted versus the 22 model-phase combinations m; 
see Sect. 4.1) and the measured R Cas effective temperatures derived from measured angular Rosscland radii i?,"^ 
(see Sect. 4.2) and JHKLM-photometry. Tabic ^ shows the link between the abscissa values m and the models and 
phases. 



4.4. Effective temperatures of R Cas 

Effective temperatures of R Cas were derived from its angular Rosseland radii i?^ (Sect. 4.2) and its bolometric flux 
using the relation 

Teff = 2341 K X (^^bol/0')'/^ (1) 

where Fboi is the apparent bolometric flux in units of 10^^ ergcm"^ s~^ and = 2i?°i is the angular Rosseland 
diameter in mas. The bolometric flux was derived from JHKLM-band observations carried out 7 days after the 
visibility observations. For cool stars such as LPVs, where most of the luminosity is emitted at near-infrared 
wavelengths, a convenient approximation for calculating bolometric magnitudes is to use a blackbody function to 
interpolate between photometric observations in the J, H, K, L and M bands. For estimating the bolometric flux we 
used the photometric JHKLM measurements of B. Yudin which were carried out with the 1.25 m telescope at the 
Crimean station of the Sternberg Astronomical Institute in Moscow seven days after our visibility observations. The 
derived bolometric flux of R Cas on 29 August 1999 was 980.6±180 xlQ-^ergcm-^ s"^ 



Table 3. Observational R Cas data and measured effective temperature. 



Star 


Date 




K 


Fbol 












[mag] 


[10"* erg/cm^ s] 


[mas] 


[K] 


R Cas 


99 August 29 


0.08 


-2.04 


980.6±180 


12.1±1.7 


2685±238 



Fig. g shows a comparison of measured and theoretical Toff^m-values. The theoretical and measured (phase 0.08) 
Toff-values are in agreement, within the la error bar, with the E models at near-maximum phases 1.0 and 1.21, 
with the P models at phases 2.0 and 4.0, and with the M models at phases 1.0 and 2.0. Table || lists the measured 
bolometric flux, the average measured angular Rosseland radius (see Sect. 4.2) and the average measured effective 



temperature of R Cas. The average measured efFeetive temperature and angular Rosseland radius is the average over 
all Tes,in- and ii^-values corresponding to m- values with phases close to the phase of our observations (i.e., m = 1, 
3, 6, 9, 11, 13, 15, 17, 19, 22 were used for averaging). 

4.5. Pulsation mode 

Adopting the above (Sect. 4.3) phase-averaged (over models M and P at near-maximum phases) linear Rosseland 
radius of the non-pulsating parent star, Rp, of 24OR0±59R0 for R Cas (derived from the measured R and the 
theoretical ratios given in Table 2), we find for the pulsation constant Q = P (Af/Af©)^/^ [R/ Rq)^'^/'^ a value of 
(5=0.116±0.042 for a 1 M© Mira with period P=430 days. This Q value agrees within the la error with the theoretical 
value ((5=0.097) for fundamental pulsation mode for IM© AGB stars with a period of -^430 days (Fox & Wood ) . 
The corresponding Q value of first overtone pulsation mode is (5=0.056. Note, however, that no direct measurement 
of a Mira mass exists and that a 20% uncertainty of M would for example result in a 10% uncertainty of Q. 



5. CONCLUSION 



We have observed the Mira star R Cas with the GI2T/REGAIN interferometer and our new infrared beam combiner. 
A uniform-disk (UD) K-band (2.1 /xm) diameter of 25.3 mas ±3.3 mas was derived from the interferograms obtained 
with baselines between 12.0 m and 13.9 m at near-maximum variability phase 0.08. The unusually high visibility 
values at baselines 18.1 m, 19.8 m, 22.0 m, and 23.8 m show that the stellar surface of R Cas is more complex than 
previously assumed. Some of the theoretical models discussed above have an unusual wing- like visibility shape (e.g., 
P76200 in HSW98). Huge convection cells (predicted by Schwarzschild^^) could also explain the observed visibility 
function. The visibility function of R Leo observed by Pcrrin et al. with the IOTA interferometer has a similar 
wing-like shape as the R Cas visibility. Perrin et al. discuss several physical phenomena which may explain the 
excess of visibility at high frequencies. 

The R Cas UD diameter measured at wavelength 1.04 /xm at near-minimum phase 0.63 is 29.9 mas ±3.0 mas 
(HBSW2000 ^^), i.e., considerably larger than the above 25.3 mas ±3.3 mas UD diameter measured at 2.1 fim at 
near-maximum phase 0.08. The reason for this difference is not known and unexpected from the theoretical models 
discussed in the previous sections. The different variability phase of the two observations and/or other unknown 
effects, for example time- variable surface structure (e.g., supergranulation^^), are possibly the cause of the diameter 
difference. 

An angular R Cas Rosseland radius i?° of 12.1 mas ±1.7mas at phase 0.08 was derived from the 2.1 /um visibilities 
measured with projected baselines between 12.0 m and 13.9 m by fitting theoretical center-to-limb variation profiles 
of five recent Mira star models (BSW96, HSW98). From the above mentioned 1.04 /xm observation an angular 
Rosseland radius R° of 16.5mas±1.7mas at phase 0.63 (HBSW2000), which is larger than the the near-maximum 
Rosseland radius reported in this paper, was obtained (HBSW2000; see above discussion). 

The effective temperature of 2685 K ±238 K at near-maximum phase 0.08 was derived from the measured angular 
Rosseland radius and JHKLM-photometry carried out only seven days after the visibility observations. 

For R Cas a good HIPPARCOS parallax (9.36 mas ±1.10 mas) is available and it is therefore possible to compare 
measured linear Rosseland and stellar filter radii with the theoretical radii of the BSW96 and HSW98 models. 
The measured radii were derived by fitting theoretical (BSW96, HSW98) center-to-limb intensity variations to the 
low-frequency visibility data. In the following table we compare measured and theoretical values: 



Measured linear Rosseland R Cas radii R: 
Theoretical linear Rosseland radii R: 
Measured linear stellar K-band R Cas radii Rm: 
Theoretical linear stellar K-band radii Ry^: 
Measured effective R Cas temperature: 
Theoretical effective temperature: 



286±51R0 
315 R© 
277±5OR0 

304 Rq 
2621±217K 



(M model*) 
(M model*) 
(M model*) 

(M model*) 
(M model*) 
(M model*) 



264±82R0 
263 R© 
252±68R0 

259 Rq 
2797±275K 
3030 K 



(P model**) 
(P model**) 
(P model**) 
(P model**) 
(P model**) 
(P model**) 



2700 K 

(* average over phases 1.0, 2.0; ** average over phases 1.0, 2.0, 3.0, 4.0; derived from visibilities measured at baselines 
between 12.0 m and 13.9 m) 

The comparison of these K-band (2.1 ^m) observations with theoretical models suggests that R Cas is well represented 
by the fundamental mode M and P model (BSW96, IISW98), whereas the above mentioned 1.04 /xm observations 



suggested first-overtone pulsation. However, observations in more filters than just one continuum, filter may be 
necessary for safely distinguishing a well-fitting model from an accidental match (cf. Ref. 16). 
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